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Abstract. In this paper, we construct a Malliavin derivative for func- 
tionals of square-integrable Levy processes and derive a Clark-Ocone for- 
mula. The Malliavin derivative is defined via chaos expansions involving 
stochastic integrals with respect to Brownian motion and Poisson ran- 
dom measure. As an illustration, we compute the explicit martingale 
representation for the maximum of a Levy process. 



1. Introduction 

If W = {Wt)t£[o,T] is a Brownian motion, then the Wiener-Ito chaos ex- 
pansion of a square-integrable Brownian functional F is given by 

T rT 

fnitl,...,tn)Widh)...W{dtn), 



(1) F = E[F] + Y^ [ ... / 



where (/n)n>i is a sequence of deterministic functions. This chaotic repre- 
sentation can be obtained by iterating Ito's representation theorem and can 
then be used to define the classical Malliavin derivative in the following way: 
if the chaos expansion of F satisfies an integr ability condition, then F is 
Malliavin-differentiable and its Malliavin derivative DF is given by 

(2) DtF = h{t) 

+ V(n + 1) / .../ fn+l{ti,...,tn,t)W{dti)...W{dtn), 
n>l 

for t € [0,r]. This derivative operator is equal to a weak derivative on the 
Wiener space; the close connection between Hermite polynomials and Brow- 
nian motion is at the hearth of that equivalence. See for instance Nualart 

m. 
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Quite recently, L0kka [H] developed similar results for a square-integrable 
pure-jump Levy process L = (-Z>t)tG[o,T] given by 



where /i — vr is the compensated Poisson random measure associated with 
L. In this setup, by mimicking the steps of the Wiener-Ito expansion, L0kka 
obtained a chaos representation property for the pure-jump Levy process L 
just as in Equation ([T]) and then defined the corresponding Malliavin de- 
rivative as in Equation Later on, Benth et al. [4] introduced chaos 
expansions and a Malliavin derivative for more general Levy processes, i.e. 
Levy processes with a Brownian component. However, in the latter, neither 
proofs nor connections with the classical definitions was given. 

Our first goal is to provide a detailed construction of a chaotic Malliavin 
derivative leading to a Clark-Ocone formula for Levy processes. We extend 
the definitions of the Malliavin derivatives for Brownian motion and pure- 
jump Levy processes to general square-integrable Levy processes. Secondly, 
we derive additional results that are useful for computational purposes. 

Our approach follows more or less the same steps as those leading to the 
Wiener-Ito chaos expansion and the chaotic Brownian Malliavin derivative, 
just as L0kka [H] did for pure-jump Levy processes. The definition of the 
directional Malliavin derivatives is different from those of Benth et al. [4]. 
The main idea is to obtain a chaotic representation property (CRP) by it- 
erating a well-chosen martingale representation property (MRP) and then 
defining directional Malliavin derivatives in the spirit of Ma et al. [16]. How- 
ever, in the context of a general square-integrable Levy process, one has to 
deal with two integrators and therefore must be careful with the choice of 
derivative operators in order to extend the classical definitions. This choice 
will be made with the so-called commutativity relationships in mind and 
following Leon et al. [13]. In the Brownian motion setup, the commutativity 
relationship between Malliavin derivative and Skorohod integral is given by 



when u is an adapted process. See Theorem 4.2 in Nualart and Vives [20] 

for the corresponding formula in the Poisson process setup. 

We will get the MRP using a denseness argument involving Doleans-Dade 
exponentials. Our path toward the CRP is different from that of Ito [9] and 
Kunita and Watanabe [12] who used random measures; see also the recent 
formulation of that approach given by Kunita [H] and Sole et al. p5]. It is 
known that the CRP usually implies the MRP and that in general a Levy 
process does not possess the MRP nor a predictable representation property. 
However, we show that the CRP and our well-chosen MRP are equivalent 
for square-integrable Levy processes. Finally, just as in the Brownian and 
pure-jump Levy setups, a Malliavin derivative and a Clark-Ocone formula 




(3) 
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are derived. As an illustration, we compute the explicit martingale repre- 
sentation for the maximum of a Levy process. 

This approach to Malliavin calculus for Levy processes is different from 
the very interesting contributions of Nualart and Schoutens [19], Leon et al. 
[T3] and Davis and Johansson [6]. They developed in sequence a Malliavin 
calculus for Levy processes using different chaotic decompositions based on 
orthogonal polynomials. Their construction also relies on the fact that all 
the moments of their Levy process exist. Many other chaos decompositions 
related to Levy processes have been considered through the years: see for 
example the papers of Dermoune [7], Nualart and Vives [20], Aase et al. [1] 
and Lytvynov [15j. 

On the other hand, Kulik [10] developed a Malliavin calculus for Levy 
processes in order to study the absolute continuity of solutions of stochastic 
differential equations with jumps, while Bally et al. [2] established an inte- 
gration by parts formula in order to give numerical algorithms for sensitivity 
computations in a model driven by a Levy process; see also Bavouzet-Morel 
and Messaoud 0. Finally, in a very interesting paper. Sole et al. ^5] con- 
structed a Malliavin calculus for Levy processes through a suitable canonical 
space. While finishing this paper, the work of Petrou [2T] was brought to 
our attention. In that paper, the same methodology is applied to obtain a 
Malliavin derivative and a Clark-Ocone formula, but the focus is on finan- 
cial applications. Since one of our goal is to give a thorough treatment of a 
Malliavin calculus for square-integrable Levy processes, and considering the 
major differences between the two papers, we think that each paper has his 
own interest. 

The rest of the paper is organized as follows. In Section 2, preliminary 
results on Levy processes are recalled. In Section 3 and 4, martingale and 
chaotic representations are successively obtained. Then, in Section 5, the 
corresponding Malliavin derivative is constructed in order to get a Clark- 
Ocone formula. Finally, in Section 6, we apply this Clark-Ocone formula to 
compute the martingale representation of the maximum of a Levy process. 



2. Preliminary results on Levy processes 

Let T be a strictly positive real number and let X = {Xt)t(z\[)^T] be a 
Levy process defined on a probability space P), i.e. X is a process 

with independent and stationary increments, is continuous in probability and 
starts from almost surely. We assume that X is the cddldg modification 
and that the probability space is equipped with the completed filtration 
(^t)tg[OT] generated by X. We also assume that the cr-field T is equal to 

This filtration satisfies les conditions hahituelles and, for any fixed time t, 
J^t- = ^t- Consequently, the filtration is continuous. This fact is crucial in 
the statement of our Clark-Ocone formula. 
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The reader not familiar with Levy processes is invited to have a look at 
the books of Schoutens [23], Protter [22] and Bertoin [5]. 

Prom the Levy-Ito decomposition (see [22], Theorem 42), we know that 
X can be expressed as 

(4) Xt = at + aWt+ [ [ zN{ds,dz)+f [ zN{ds,dz) 

Jo J\z\>l Jo J\z\<l 

where a is a real number, a is a strictly positive real number, is a stan- 
dard Brownian motion and N is the compensated Poisson random measure 
associated with the Poisson random measure N. The Poisson random mea- 
sure is independent of the Brownian motion W. Its compensator measure 
is denoted by A x i/, where A is Lebesgue measure on [0, T] and u is the Levy 
measure of X, i.e. u is & a-finite measure on M such that i^({0}) = and 

/ {lAz'^)iy{dz) < oo. 

Therefore the compensated random measure is defined by 

iV([0, t]xA) = N{[0, t]xA)- tiy{A). 

This measure is equal to the measure — vr mentioned in the introduction. 

Finally, let V be the predictable cr-field on [0, T] x n and B{R) the Borel 
cj-field on R. We recall that a process ip{t,z,uj) is Borel predictable if it is 
(V X ;B(R))-measurable. 



2.1. Square-integrable Levy processes. When the Levy process X is 
square-integrable, it can also be expressed as 



(5) 



Xt = fit + aWt + 







z N(ds, dz), 



where ^ = ]E[Xi]. Indeed, in Equation (|4| we have that 



a = E 



zN[dt,dz) 



J|^|>1 



so ]E[Xf] is finite if and only if 



z'^v{dz) = E 



-'|2|>l 



zN{dt,dz) 



is finite. Note that in general ijl ^ a. 
Here is a consequence of Ito's formula. 
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Lemma 2.1. If h belongs to L'^{[0,T],X) and if {t,z) ^ e^^*'^) - 1 belongs 
to L\[0,T] X M, A X ly), define Z = {Zt)te[o,T] by 

(6) Zt = exp\ [ h{s)W{ds)-l f h^{s)ds+ [ [ g{s,z) N{ds,dz) 
[Jo 2 7o Jo Jr 

- f j (e»(^'^) - l) v{dz)ds\ . 

The process Z is a square-integrable martingale if and only if — 1 is an 
element of L'^{[0,T] x M,A x z/). 

Proof. Prom the assumptions, we have that g belongs to -^^^([O, T] x M, A x z^) 
and that Z is a well-defined positive local martingale. Then, if E[Zt] = 1, 
it is a martingale. From Ito's formula, we also have that Z is the solution of 

dZt = Zt- h{t) W{dt) + Zt- [ (es(*'^) - 1) N{dt, dz), Zq = 1. 

Let {Tn)n>i be the fundamental sequence of stopping times of Z. Since W 
and N are independent. 



nzirj = i+^' ' 



Z^h^{s) ds 





-E 



/ Z^ / (e»(^'^) - if u{dz)ds , 
Jo Jr J 

for every n > 1. Taking the limit when n goes to infinity yields 

(7) E[Zf] = 1 + /* h'^{s) ds + /* E[Z^] [ (e^(^'^) - 1)^ u{dz)ds. 

Jo Jo Jr 

If we define G{t) = h?{t) + ^(e^'^*''') - 1)^ v{dz), then the function t ^ E[Z|] 
is the solution of 

F'it) = G(i)F(i), F(0) = 1. 

Hence, 

(8) E[Z^] = exp I r h'^{s) ds+ f I (es(*'^) - 1)^ u{dz)ds 

I Jo ^0 JR 

and the statement follows. □ 

For h e L2([0,r],A) and e^? - 1 G L2([0,T] x M,A x ;/), the proc_ess Z is 
the Doleans-Dade exponential of the square-integrable martingale (-/W^t)te[o,T] 
defined by 

M^ = f h{s) W{ds) + f f (e^(^'^) - 1) N{ds, dz). 
Jo Jo Jr 

In the literature, this is often denoted by Z = £{M), the stochastic expo- 
nential of M. 
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2.2. A particular choice for g. If g is an element of -^^^([0, T] x M, A x i/), 
then — 1 is not necessarily square-integrable. One way to circumvent this 
problem is to introduce the bijection 7: IR ^ (—1, 1) defined by 

(9) 7(^) = 

Note that 7 is bounded. Hence, if h is square-integrable on [0, T] and if g 
is of the form g{t,z) = g{t)^{z), where g G C([0, T]), i.e. ^ is a continuous 
function on [0, T], then Z is square-integrable by Lemma [2711 

The idea of introducing the function 7 is taken from L0kka |14| . In that 
paper, it is also proved that the process {Nt)ti^[o^T] defined by 

(10) Nt= [ [ zN{ds,dz) 



and the process {Nt)ti^[Q^T] defined by 







Nt= / -f{z)N{ds,dz) 



generate the same filtration. Since 

for every t £ [0,T] (see Lemma 3.1 in [SSj), we have the following lemma. 
Lemma 2.2. For every t £ [0,T], 

As a consequence, T = V . 

This means that the processes Xt = [it^aWt^^t and Xt = /it+aWj + A't 
both generate the filtration {Tt)t^\^Q^T\- 

3. Martingale representations 

Assumption 3.1. For the rest of the paper, we suppose that X is a square- 
integrable Levy process with a decomposition as in Equation ^ . 

In general, a Levy process does not possess the classical predictable repre- 
sentation property (PRP), i.e. an integrable random variable F (even with 
finite higher moments) can not always be expressed as 

t-T 



F = E[F]+ / UtdXt, 
Jo 



where u is a predictable process and where the stochastic integral is under- 
stood as an integral with respect to a semimartingale. However, a martingale 
representation property exists for square-integrable functionals of X. It is 
a representation with respect to W[dt) and N{dt, dz) simultaneously. This 
result can be found as far back as the paper of Ito [9j- In this section, we 
will provide a different proof. But first, here is a preparatory lemma. 
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Lemma 3.2. The linear subspace o/L^(r2,^, P) generated by 

{Y{h,g) \ heL\[0,T],X), geC{[0,T])} , 
where the random variables Y{h,g) are defined by 

(11) Y{h, g) = exp I r h{t) W{dt) + T f g{t)^{z) N{dt, dz)] 

Uo Jo Jr ) 

is dense. 



Proof. We adapt the proof of Lemma 1.1.2 in the book of Nualart [18]. Let 
X be a square-integrable random variable such that 

E[Xy(/i,5)] =0 

for every h G L2([0, T], A) and g G (7([0, T]). Let W{h) = h{t) W{dt) and 
^(9) = J^g{th{z)N{dt,dz). Hence, 



E 







Xexp<^^(aiVF(/ii) + W5,)^ 
U=l 

for any n > 1, any {oi, . . . , a„, 6i, . . . , 6„} C M and any (sufficiently inte- 
grable) functions {hi, . . . , hn, gi, ■ ■ ■ , gn}- Then, for a fixed n and fixed func- 
tions {hi, . . . , hn, gi, ■ ■ ■ , gn}, the Laplace transform of the signed measure 
on S(M") X ^(M") defined by 

{A, B)^E [xiA{W{hi), Wihn))lB{N{gi), N{gn)) 

is identically 0. Consequently, the measure on ^ = Tt defined hy E 
E [XI^;] vanishes on every rectangle ^ x i? if it is a pre-image of the M^"- 
dimensional random vector 

W{hi),...,W{hn),N{gi),...,N{gn) 

By linearity of the stochastic integrals, this is also true for random vectors 
of the form 



W{hi),...,W{hn),N{gi),..., Nigra)) , 



when m and n are different. Since is generated by those random vectors, 
the measure is identically zero and X = 0. □ 

We now state and prove a Martingale Representation Theorem with re- 
spect to the Brownian motion and the Poisson random measure simultane- 
ously. 

Theorem 3.3. Let F G L^{n,T,F). There exist a unique Borel predictable 
process G L^(A x x P) and a unique predictable process (/) G -L^(A x P) 
such that 

(12) F = E[F]+ / ct){t)W{dt)+ [ [ il}{t,z)N{dt,dz). 
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Proof. For h G L^([0,T],A) and g G C([0, T]), we know from the proof of 
Lemma 12.11 that 

yt = exp|/" h{s)W{ds)-l [ h^{s)ds+ I [ g{s)j{z) N{ds,dz) 
I .In ^ .In .In .Ins 



,9{s)l{z) _ ^ _ g{^s)^{^z)^ u{dz)ds 



is a solution of 



(13) Yt = l+ f Y,^h{s) W{ds) + f j (e^^(-hW _ i) N{ds, dz) 
Jo Jo Jr ^ ' 



on [0, T]. Hence, Yt admits a martingale representation as in Equation (fT2]) 
with (\){t) = Yt-h(t) and ^lj{t,z) = Yt- (e^^W^^^) - l). These two processes 
are predictable. Note that 

Yt = y(/i,5)e-^^(^'5) 



where 



OriKg) = \f h\t) dt+ f [ (e^W^(^) - 1 - g{tHz)) v{dz)dt. 
^ JO Jo Jr^ ' 

Since Oxih^g) is deterministic, Y{h,g) also admits a martingale representa- 
tion as in Equation (fT2]) but this time with 

(j){t) = Yt-h{t)e^^^''^'^^ and ^(t, z) = Yt- (^e^^W^^^) - l) e^^^^'^l 

Therefore, the first statement follows by a denseness argument. Indeed, from 
Lemma [321 since F is square-integrable, there exists a sequence (-Fn)n>i of 
square-integrable random variables such that F„ tends to F in the L'^{Q)- 
norm when n goes to infinity. Moreover, the Fn's are linear combinations of 
some Y(h,g)'s. Then, for each term in this sequence there exist cpn and ^jJn 
such that 

Fn=E[Fn]+ [ (t>n{t)W{dt)+ [ [ i^^^it, z) N {dt, dz). 

Jo Jo Jr 

Also, since 

rT 



E[Fn - Fmf = E 



E[Fn -Fm]+ / (Mt) - <Pmit)) W (dt) 





+ [ [ {Mt,z)-i^m{t,z))N{dt,dz) 

Jo Jr 

rT 

{E[Fn-Fm]f+ / E[cl)n{t) - (l)m{t)? dt 

Jo 

+ [ f E[llJn{t,z)-i)m{t,z)fv{dz)dt, 

Jo Jr 
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we get that {4>n)n>i and ('0n)n>i are Cauchy sequences. It follows that there 
exist predictable processes i/j ^ L'^ {X x v x F) and </> G L^(A x P) for which 
the representation of Equation (fT2]) is verified. 

We now prove the second statement. If F admits two martingale repre- 
sentations with 01 , <^2 , V'l ) ^2 (these have nothing to do with the previous 
sequences), then by Ito's isometry 

= \\<Pl - Ml^XxP) + IIV'l - ^2\\l2(^xxuxV) 

and then (/>! = 4>2 in L^(A x P) and tpi = ip2 in L^(A x x P). □ 

Remark 3.4. From now on, we will refer to this martingale representation 
property of the Levy process X as the MRP. 

4. Chaotic representations 

We now define multiple integrals with respect to W{dt) and N{dt, dz) 
simultaneously and define Levy chaos as an extension of Wiener-Ito chaos. 
Then, we show that any square-integrable Levy functional can be represented 
by a chaos expansion. We refer the reader to the lecture notes of Meyer [T7] 
for more details on multiple stochastic integrals. 

4.1. Notation. In the following, we unify the notation of the Poisson ran- 
dom measure and the Brownian motion. Thus, the superscript (1) will refer 
to Brownian motion and the superscript (2) to the Poisson random measure. 
This is also the notation in |^. 

Let X = [0, T] X M. We introduce two (projection) operators IIi : A" — > 
[0,T] and II2 . X ^ X defined by Ili{t,z) = t and Il2{t,z) = {t,z). Conse- 
quently, Hi ([0,r] X R) = [0,r] and U2 ([o,r] X M) = [o,r] X M. 

For n > 1, t € [0,T] and (zi,...,i„) € {1,2}", we also introduce the 
following notations: 

(14) S„(t) = {(ti,...,t„) G [0,r]" \ ti < ■■■ <tn<t}; 

and 

S(n,...,i„)([0,t] xM) 

= {{xi,...,xn) €Ui,{A:) X ••• X Uijx) I ni(xi) < ••• < ni(x„) < t}. 

Consequently, T,n{T) = j^-)(A') when ik = 1 for each k = 1,2,... ,n. 

If / is a function defined on x ■ ■ ■ x we write f{xi, . . . ,Xn), 

where Xk G for each /c = 1, 2, . . . , n. If r/i = A and 772 = A x i^, 

let L2 (S 

(ii,...,iri)i^)) be the space of square-integrable functions defined on 
S(ii,...,i„)('^) and equipped with the product measure r/j^ x • • • x r]i^ defined 
on Ui[iX) X ••• X HijA"). 
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4.2. Multiple integrals and Levy chaos. Fix n > 1 and G 
{1, 2}". We define the iterated integral J(ji,...,i„)(/), for / in [T,^i^^_^i^){X)) , 
by 

J{ii,...,in)if) 

= [ ...I f{xi,...,Xn)M^''\dxi)...M^"^\dXn) 

Jlli„{[0,T]xR) Jlli^{[0,t2~]xM.) 

where M^\dx) equals W{dt) if j = 1 and equals N{dt, dz) if j = 2. The i\ in 

J{ii,...,in) stands for the innermost stochastic integral and the in stands for the 
outermost stochastic integral. For example, if = 3 and (ii, ^2, is) = (1, 1, 2), 
then 



•^(1,1,2) (/) 



T 



fitl^t2,it3,Z3))W{dh))W{dt2) 



N{dt3,dZ3). 



As n runs through N and runs through {1,2}'^, the iterated 

integrals generate orthogonal spaces in that we would like to call Levy 

chaos. Indeed, since 



ii([o,r]x 

and 



/ f{x)M^\dx) 



! g{x)M^^\dx) 

in^-([0,T]xR) 



are independent \i i ^ j and both have mean zero, using Ito's isometry 
iteratively, we get the following proposition. 

Proposition 4.1. // / G L2(E(,,,...,,„)(A')) and if g & L2(E(,-,, (A")), 
then 

E[J(n,...,i„)(/)Joi,...j„)(5)] 

^ f(/,5)L2(S(i^,...,i„)(Ar)) if{ii,--.,in) = (ii,---,Jm); 
1 if not. 

We end this subsection with a definition. 

Definition 4.2. For n > 1 and (ii, . . . , in) G {1, 2}", the {ii, . . . , in)-tensor 
product of a function h defined on [0, T] with a function g defined on [0, T] xM 
is a function on Hj^ (X) x • • • x Hj^ (X) defi,ned by 



{h(S){i^,...,i„)g) {xi,...,xn) = Y[ h{Ui{xk))'^ g{U2{xk)y'' V 

l<fe<n 

For example, 

(^<^(i,i) g) {s,t) = h{s)h{t) 
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is a function defined on [0, T] x [0, T] and 

[h ®(i,2,i) g) {r, {s, y),t) = h{r)h{t)g{s, y) 
is a function defined on [0,T] x ([0,T] x M) x [0,T]. 

4.3. Chaotic representation property. For the rest of the paper, we will 
assume that E{ii,...,i„) means E(»i,...,i„)e{i,2}" ■ 

Recall that Z = (^t)te[o,T] was defined in Equation ([6|) by 

Zt = expl [ h{s)W{ds)-- [ h'^{s)ds+ [ [ g{s, z) N{ds,dz) 
[Jo 2 Jq Jo Jr 

- u{dz)di 

Lemma 4.3. Let h G L'^{[0,T]) and - 1 £ L'^{[0,T] x R, X x u) . Then, 
Zt admits the following chaotic representation: 

oo 

(15) Zt = 1 + Y^ '^(n,...,.„)(^ (e^ -I))- 

n=l 

Proof. We know from the proof of Lemma 12.11 that Zt is square-integrable 
and that 



(16) 



Zt = 1+ [ Zt^h{t)W{dt) + [ [ Zt_(e^^(*'^) - l)N{dt,dz). 
Jo Jo Jr 



Let (l)^^\t) = Zt-h{t) and (p^'^\t,z) = Zt_(e^'(*'^) - 1). We now iterate 
Equation ([T6]) . Consequently, 



zt=i + r f^^\t)w{dt)+ r [ f^^Ht,z)Nidt,dz) 

Jo Jo Jr 

[ Zs-h{s)h{t)W{ds)W{dt) 
'0 Jo 

+ [ [ [ Zs-{e3^''y^ -l)h{t)N{ds,dy)W{dt) 

Jo Jo Jr 



T r rt- 



JR Jo 
T r ft- 



JR Jo 



Z,_/i(s)(e^?(*'^) - 1) W{ds) N{dt, dz) 

Z^^{e3^'^y) - l)(es(*'") - l)N{ds,dy)N{dt,dz) 
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where f^^\t) = h{t) = (/i (e^ - 1)) (t) and f^^\t,z) = e3(*'^) - 1 = 
{h 0(2) (e^ - 1)) (t, z). Then, after n iterations, we get 

n-l 
k=l 

(n,...,i„) -^n.^ {[o,T] xR) Jn,^ {[o.ta-] xR) 

where = /j (e^ - 1) and where = _^_(/^ ..,«„) 

(e^ — 1)). This means that we can define a sequence {ipn)n>2 in L'^{Q) by 

(ii,...,i„)-^n,„([0,T]xR) Jn,j([0,t2-]xR) 

M(*i)((ixi)...M(^")(d3;„). 



Prom Proposition 14.11 

n-l 

k=l («!,.. 

for each n > 2. Hence we get that 

oo 

is a square-integrable series and that there exists a square-integrable random 
variable ijj such that il)n tends to if) in the L^(ri)-norm. Consequently, it is 
enough to show that ■0 = 0. Since /(^i----'*") = h (8)(ji,...,i„) (e^ — 1), using 
Proposition 14.11 once again, we get that 



{il,...,in) 



E E 11^ (e^ - i)iii2(E(,^,...,„)W)' 

fc=0 



i|=fc 



where |i| = . . . = X^"=i(2 — ij) stands for the number of times the 
function h appears in the tensor product. Note that when \i\ = k there 
are (^) terms in the innermost summation. Since (e^ — 1)^ is a 

(ii, . . . , i„)-tensor product, the function given by 

E (e'-i) 
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E 



fc=0 ^(H>--->>n)('^) '- 

\i\=k 



%i • • • drjin 



1 " 

A:=0 
1 " 



lil=fc 



. . . dr]i^ 



1 



A:=0 



n 



\2k 



\p9 _ 1 ||2{«-fc) 
^ ^IIl2(Ax!^) 



li2(A) + W^^ ~ ^Wl^iXxu) 



Prom Equation ([8|), we know that 

E[Z|] = exp{||/i||i,(,) + ||e^-l||i,(,,,)}. 
This means that -0 = and the statement follows. 



□ 



We are now ready to state and prove the chaotic representation property 
of the Levy process X. The previous lemma and the idea of its proof will be 
of great use. 

Theorem 4.4. Let F G L'^{n,J^,F). There exists a unique sequence 
whose elements are respectively in I? (S^j^^ ^^■((A:')), such that 



(17) 



Consequently, 



n=l («!,.. 



(18) E[F']=E'[F] + J2 E 



f{ii,...,in) ||2 



=1 



Proof. Prom Theorem 13.31 we know there exist a predictable process (p^^^ G 
L^(A X P) and a Borel predictable process 4>^'^^ G L^(A x x P) such that 



F = E[F]+ / </.W(t)VF(dt) + / / ^'^^\t,z)N{dt,dz). 



T 



Using Ito's isometry, it is clear that 



(1)||2, 



lL2(AxI 



+ 



/,(2)||2 



L2(Axi/xIf 
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For almost all t G [0,r], G L2(J1, jf^, P) and then from Theorem EJ 

there exist processes tp'^'^'^^ and i?!)^^'^) such that 



<A«(f)=E[</.«(t)]+ / <l^^^^^\t,s)W{ds)+ / / (t>^^^^\t,s,y)N{ds,dy). 
Jo Jo Jr 

Similarly, for almost all {t,z) G [0,T] x M, (p'-^\t,z) G L2(17,JS,P) and 
(/>(2) {t, z) = E[,/.(2) (t, z)] + r ,^(2,1) (t^ s) VF(ds) 

+ /" / ^^^^^\t,z,s,y)N{ds,dy). 
Jo Jr 

Consequently, 

F = E[F]+ [ g^^\t)Widt)+ [ [ g^^\t, z) N{dt,dz) 
Jo Jo Jr 

+ [ [ (l)^^'^Ht,s)W{ds)W{dt) 

Jo Jo 

+ f f f (p^^'^Ht,s,y)N{ds,dy)W{dt) 

Jo Jo Jr 

+ [ If cl)'^^''^\t,z,s)W{ds)N{dt,dz) 
Jo Jr Jo 

+ [ [ [ I (l)^^'^\t,z,s,y)N{ds,dy)N{dt,dz). 
Jo Jr Jo Jr 

where g^^\t) = E[(j)^^\t)] and g^'^\t,z) = E[4>^'^\t, z)]. After n steps of this 
procedure, i.e. after n iterations of Theorem 13.31 we get as in the proof of 
Lemma 14.31 that 

n-l 
k=l 

where g ... , for each 1 < /c < n-l and (ii, . . . G 

{1, 2}^, where 

i^n= J2 f - I </.(*^'-'*")(xi,...,x„) 

-^n^n ([0,T] xM) JU,^ ([0,t2H xR) 

M^'^\dxi)...M^'-\dxn), 

and where g 1,2 (^^^^ x • • • x ryj„ x P), for each {ii, ■ ■ ■ ,in) G {1)2}". 

Prom Proposition 14.11 

n-l 

fc=l (ii,...,ife) 
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for each n > 2 and 

oo 

E E ^(n,...,.0(/^"'-'^"^) 

n=l 

is a square-integrable series. Consequently, we know that there exists a 
square-integrable random variable ijj such that ipn tends to ip in the L^(r2)- 
norm. It is enough to show that ip = 0. Using the argument leading to 
Proposition [4]T1 i.e. the fact that two iterated stochastic integrals of different 
order are orthogonal, we get that for a fixed n > 2, 



L2(n) 



for every l<k< n-1, (i^, . . . , 4) g {1, 2}'^ and /(^i.-.^fe) ^ L\E^i^_,^){X)). 
Thus, 



(19) (j(n,...,.„)(/(*^'-'^")),V') 



= 

L2{n) 



for every n>l, (ii, . . . , i„) G {1, 2}" and /(^^■•••'*") G L2(S(,^_,^)(;f)). 

We now assume that g = where g belongs to C([0, T]). Using Equa- 
tion (fT9]l . we have that ip is orthogonal to each random variable Y{h,g) 
defined in Equation ([TT]) since from Lemma 14.41 they each possess a chaos 
decomposition. We also know from Lemma [32] that these random variables 
are dense in L'^{^, so ip = 0. This means that every square-integrable 

Levy functional can be express as a series of iterated integrals. The state- 
ment follows. □ 

Remark 4.5. From now on, we will refer to the chaotic representation prop- 
erty of Theorem \4-4\ ^5 the GRP. 



Remark 4.6. As mentioned before, in general the CRP implies the MRP. 
Indeed, if F is a square-integrable Levy functional with chaos decomposition 

00 

n=l («!,.. 

then ^ ^ 

F = E[F]+ / <l){t)W{dt)+ [ [ ^p{t,z)N{dt,dz), 
Jo Jo Jr 

with 

00 

n=l {ii,...,i„) 
00 

i;{t,z) = f^'\t,z) + J2 E (/^'^'•••''"''H-,(i,^))IlE„W 

n=l (ii,...,i„) 

This last remark, together with our journey from the MRP of Theorem l3.3l 
to the CRP of Theorem 14.41 yields the following interesting proposition. 
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Proposition 4.7. For a square-integrable Levy process, the MRP and the 

CRP are equivalent. 

4.4. Explicit chaos representation. In the next proposition, we compute 
the exphcit chaos representation of a smooth Levy functional. 

Proposition 4.8. Let f he a smooth function with compact support in M'^, 
i.e. let f G C^i'R''), and let tj belong to [0, T] for each j = l,...,k. Then, 

oo 

f{Xt„...,Xt^)=E[f{Xt„...,Xa + Yl E %,.....)(/^^"-""^)' 

n=l (ii,...,i„) 

where 

f('^'-''-\ni,{si,Wi),...,Ui„{Sn,Wn)) 

■^^^ l<j<n 

with 

4>{x) dx = ¥{Xt G dx}, 
where Xt = [Xt^ , • • • , Xt^), and with 

e^'^ = yii[o,ti](s) + --- + yfci[o,tfc](«)' 

fort = {ti,...,tk) andy = {yi,...,yk)- 

Proof. We follow an idea in [Hj and use Fourier transforms. By the Fourier 
inversion formula 

(20) fix) = (2^)-'=/2 [ f{y)e'^^^y^ dy 



where / is the Fourier transform of / and {x, y) denotes the scalar product 
in M'^' of a; = (xi, . . . , x^) and y = (yi, . . . , y^). Let (i)~{x) = (j){—x). If we 
define F{x) = E [f{Xt^ +Xl,...,Xt^^ + Xfc)], then 

F{x) = -{f * <P'){x) 

and also 

F{x) = -{2^r^/^ [ f{y)^-{y)e'^-'y'> dy. 



Therefore, we have the following equality: 

nf{Xt,,...,Xt,)] = -{2T,)-^'^ [ f{y)P{y)dy. 
From Equation (f20l) and Equation ([5]), we have that 

fiXt,,. . . ,XtJ = (2vr)-^'/2 / /(y)e^<^*'J^) dy 

= {2tt)'''/^[ f{y)e'f''^^'y^Y^'y dy 
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where 

y*'^ = exp I /" iail^y W{ds) + [ [ iz^l^^ N{ds, dz)\ . 

[Jo Jo Jr J 

Hence, 

(21) f{Xt„...,Xt,) = {27ry'/' [ f{y)e^^^''y^Z''yE[Y''y]dy 
where 

Z''y = expj r ia^i'y Wids) + ds 
[Jo ^ Jo 

+ r [ iz^i'y Nids,dz) - [ [ {e''^''' -l)u{dz)ds\. 



From Lemma [4.31 we know that 

oo 

^*'' = 1 + E E ^(n,....„) ((^^e*-^) (e"^*" - 1)) • 

n=l 

On the other hand, 

= -e-'^'^^'y^^i-y) 

Then, using Equation (f2T]l and by Lebesgue's dominated convergence theo- 
rem, 

= -(27r)-'=/2 / f{y)^{-y)dy 



(27r)-'=/2 / f{y)^{-y) 



oo 

X E E ^(n,.....) ((^^^'') ^(^1,.....) (e^^^"' - 1)) dy 

n=l («!,.. .,j„) 
oo 

= E[/(X,„...,X,J]+J] ^ J(n,....„)(/^^'-'*"^ 

n=l («i,...,i„) 

where 

/(^i'-'^")(ni,(si,?x;i),...,ni„(s„,u;„)) 

= -(2vr)-'^/2 / f{y)^{-y) ((iae*'^) (e"?"^ - 1)) dy. 

The statement follows from Definition 14. 2[ □ 
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5. Malliavin derivatives and Clark-Ocone formula 

Before defining the Malliavin derivatives, we introduce a last notation: for 
n > 1 and 1 < k < n + 1, define 

= {{h,...,tn) G [0,r]" I ti < ••• < tk-l <t<tk<--- <tn}, 

i.e. t is at the k-th position between the tj's, where to = and tn+i = T. 
Note that = S„(t), where the latter was defined earlier in Equa- 

tion (fT4l) . In a multi-index (ii, . . . , we will use ik to denote the omission 
of the k-th. index. 

We want to define two directional derivative operators in the spirit of 
Leon et al. [13| : one in the direction of the Brownian motion and one in 
the direction of the Poisson random measure. If F = J{i^^,,,^i„){f), then we 

would like to define D^^^ F and D^^^F as follows: 

n 

^=1 k-l n-k 

and 

n 

k=l k~l n—k 

where J(7)(/) = /. 

Definition 5.1. Let D^'^ = D^^) n D^^)^ ^^here if j = I or if j = 2, D^-'') is 
the subset of L'^{Q,J^,F) consisting of the random variables F with chaotic 
representation 

oo 
n=l 

such that 

oo n „ 

n=i{ji,...,i„)fc=i -^n^W 

where the inside norm is the Lp'iTi,. r~ . AX))-norm. 

Prom Theorem 14.41 it is clear that ©^'^ is dense in L^(r2), since every 
random variable with a chaos representation given by a finite sum belongs 
to 
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Definition 5.2. The Malhavm derivatives D^^) : ]D)(i) ^ ([0,r] x and 
£)(2) . i])(2) ^ ([0,T] X M X J^) are de/?ned % 

cx) ri 

n=l (ii,. ..,?„) fc=l n-fc 

anc? 

oo n 

+ E E E%.=2}^(n,...,i.,....„)(/^^^'-'"^(^'(M),^)iis._^(,)) 

n=l (ii,...,i„) fe=l fc-l n-fc 

i/ F = E[F] + E~ 1 E(n,...,^„) -^(n.-.-in) (/(^^■•••'^")) is in D(i) or ©(2) . 

Remark 5.3. For an iterated integral, the Malliavin derivatives have a 
property similar to the classical commutativity relationship. Indeed, if F = 
J{h,...,in)if)> then 

if in = ^ and 

4'> = '^{n,...,^.-i)(/(-'(*'^))Il^n-i(t)) 

■Jt JR 

if in = 2. A similar result holds for D^^^F. 

Remark 5.4. If F = ^[F]+Yl'i^=i Jn{fn), where Jn = </(i,...,i) is the iterated 
Brownian stochastic integral of order n, then 

oo n 

n=2 k=l 

oo 

= fl{t) + Y,Jn-lifn{;t)), 

n=2 

because Y1^=i^t,'' j(t) ~ ^[o,T]{t)- This is the classical Brownian Malliavin 
derivative of F. The same extension clearly holds for the pure-jump case if 
the 1 's are replaced by 2 's. 

The definitions of B'^^) and D^^^ come from the fact that we want the 
codomains of D^^^ and ^e ([0,r] x O) and ([0,r] x M x J^) 

respectively. We finally define a norm for DF = {D(^^F,D'-'^^F) in the 
following way: 

\\DFf = p«F||i2(;,,p) + p(')i^||i2(Ax.xP)- 
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This is a norm on the product space L^(A x P) x L^(A x x P). 

5.1. Properties and interpretation of the Malliavin derivatives. We 

begin this section with a result concerned with the continuity of D. It is an 
extension of Lemma 1.2.3 in Nualart [T^. The proof is given in Appendix lAl 

Lemma 5.5. If F belongs to L'^{Q), if {Fk)k>i is a sequence of elements in 
D^'^ converging to F in the L'^{Vl)-norm and i/sup^>;^ ll-^-^fcll < oo, then F 
belongs to B^'^ and {DFk)k>i converges weakly to DF in L^(A x P) x L^(A x 
V X P). 

There is a similar and stronger result stated in [14] (Lemma 6); however 
we are unable to fill a gap in its proof. 

The choice for the definitions of the Malliavin derivative operators was 
made to extend the classical Brownian Malliavin derivative as well as the 
Poisson random measure Malliavin derivative in a wider sense than Re- 
mark[531 As mentioned in the introduction, the classical Brownian Malliavin 
derivative can be defined by chaos expansions and as a weak derivative. In 
Nualart and Vives [20], it is proven that for the Poisson process there is an 
equivalence between the Malliavin derivative defined with chaos decomposi- 
tions and another one defined by adding a mass with a translation operator. 
This last result was extended by L0kka [14] to Poisson random measures. 
But now we will follow an idea of Leon et al. [13] to prove that our deriv- 
ative operators are extensions of the classical ones. Their method relies on 
the commutativity relationships between stochastic derivatives and stochas- 
tic integrals and on quadratic covariation for semimartingales; consequently, 
it is easily adaptable to our more general context. The details are given in 
Appendix [0 

Theorem 5.6. OnD^^^ the operator D^^^ coincides with the Brownian Malli- 
avin derivative and on D*-^^ the operator D^"^^ coincides with the Poisson 
random measure Malliavin derivative. 

Hence, if F S B*^^-', all the results about the classical Brownian Malliavin 
derivative, such as the chain rule for Lipschitz functions, can be applied to 
see Nualart [18] for details. But this is also true for the Poisson 
random measure Malliavin derivative. For example, an important result in 
L0kka [H] is that if F = g{Xt^ ,...,XtJ £ B^^) and 

{t,z)^ g{Xt,+zIio^t^]it),...,Xt„+zIio^tr.]it)) -g{Xt„...,XtJ 
belongs to L^(A x x P), then 

F = g {Xt, + zI[o,i,] {t),...,Xt„+ zI[o,i„] {t))-g{Xt,,...,XtJ. 

This is the adding a mass formula. Consequently, it also applies in the 
context of a square-integrable Levy process. 
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5.2. A Clark-Ocone formula. We now state and prove a Clark-Ocone 
type formula. This formula gives explicitly the integrands in the martingale 
representation of Theorem 13. 3l for a Malliavin-differentiable Levy functional. 
It is interesting to note that no particular property of the directional deriva- 
tives are needed. 

Theorem 5.7. If F belongs to B^'^, then 

F = E[F] + [ E[Df [ Tt] W{dt) + [ [ E[DfjF [ J^t] N{dt, dz). 
Jo Jo Jr 

Proof. Suppose that F has a chaos expansion given by 

oo 

F = E[F]+Y^ 

n=l 

If for example we consider the derivative operator Z)^^^ , then from Remark l4.6l 
we have to show that 

(22) E[Z)gF|^i] 

oo 

= f('\t,z)+Y, E ^(n,....„)(/^"'-'^"''^(-,(i,^))lE„w)- 

n=l 

If ik = 2, then 

_ Jo if A; = l,2,...,n-l; 

because when k = 1,2, ...,n — 1 the outermost stochastic integral in the 

(2) 

iterated integral X. 7 ■ n starts after time t. By the definition of D) !F , 
this implies that Equation ([22]) is satisfied. The same argument works for 
the derivative operator D^^^ and thus the result follows. □ 

6. Martingale representation of the maximum 

Our main goal was to provide a detailed construction of a chaotic Malliavin 
derivative and a Clark-Ocone formula. Now, to illustrate the results, we 
compute the explicit martingale representation of the maximum of the Levy 
process X. 

For < s < t < T, define Ms,t = s\\^s<r<tXr and Mt = Mo^f If 
E[Mt'] < 00, then one can show that 

POO 

(23) E[Mt \J^t] = Mt+ / FT-t{z) dz, 

JMt-Xt 

where Fs{z) = F{Ms > z}; see Shiryaev and Yor |24] and Graversen et al. 
[8]. We will use this equality to prove the next proposition. 
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Proposition 6.1. If X a square-integrable Levy process with Levy-Ito de- 
composition 

Xt = fit + aWt + [ [ z N{ds, dz), 
Jo Jr 

then its running maximum admits the following martingale representation: 

Mt = E[Mt] + [ (pit) W{dt) + [ [ z) N{dt, dz) 
Jo Jo Jr 

with 4>{t) = aFT-t{a) and i){t, z) =E [{Mr-t + z- a)+] - /~ FT-t{x) dx, 
where a = Mt — Xt. 

Proof. Since X is a square-integrable martingale with drift, from Doob's 
maximal inequality we have that Mt is a square-integrable random variable; 
see Theorem 20 in Protter [22]. Let {tk)k>i be a dense subset of [0,T], let 
F = Mt and, for each n > 1, define F„ = maxjXt^ , . . . , } . Clearly, 
{Fn)n>i is an increasing sequence bounded by F . Hence, F„ converges to F 
in the L^(r2)-norm when n goes to infinity. 

We want to prove that each F„ is Malliavin differentiable, i.e. that each 
Fn belongs to B^'^ = D^^) nB^^). This follows from the following two facts. 
First, since 

(xi, ... ,Xn)^ max{a;i, . . . ,x„} 

is a Lipschitz function on M" and since D^^^ behaves like the classical Brow- 
nian Malliavin derivative on the Brownian part of -F„, we have that 

n n 
k=l k=l 

where Ai = {F„ = Xt,} and Ak = {Fn + Xt„...,Fn / Xt^_„Fn = 
Xt,,} for 2 < A; < n. This implies that sup„>x \\D^^^ Fn\\i2i^[Q^T]xn) ^ cr^T. 
Secondly, since D^'^^ behaves like the Poisson random measure Malliavin 
derivative on the Poisson part of F„ , we have that 

<|I)Sf„| =1 max{Xt, + z^t<t^}, • • • ,^t„ + %<t„}} - i^n| < \z\, 
where the equality is justified by the following inequality: 

||max{Xt, +zI{t<t^},...,Xt„ +z%<t„|} -F„||^2([o,T]xRxn) 

<T I z^u{dz). 
Jr 

Indeed, if z > 0, then 

< max {Xt, + zl{t<ti}, ■■■,Xt„+ zl{t<t„}} - Fn < z, 
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and, if 2; < 0, then 

= Fn + mm{-Xt-^ + \z\I{t^t^y,...,-Xt„ + |z|I{t<t^}} 
= min{F„ - Xt^ + \z\I{t^t^}, . . . ,F„ - Xt„ + [z|I{t<t^}} 



< z 



This imphes that sup„>i \\D^^'^ Fn\\Lmo,T]xRxn) < T f^z"^ u{dz). 

Consequently, sup„>i |lDF„|p < T(c7^ + z^ u{dz)) and by Theorem 
we have that F is Malhavin differentiable. By the uniqueness of a weak hmit, 
this means that taking the hmit of D^^^ Fn when n goes to infinity yields 

Dl'^F = aI[o,r]{t), 

where r is the first random time when the Levy process X (not the Brownian 
motion W) reaches its supremum on [0,T], and 

Z)gF= sup (X, + zl{t<,}) - Mt. 



o<s<r 



Hence, 



E 



a^{Mt < Mt,T I ^i} 
a¥{MT-t > a}, 



where a = Mt — X^. Since Mt^x — Xt is independent of J^t and has the same 
law as Mt-i, then using Equation ([23]l we get that 



E 



E 



sup [Xs + zl{t<s}) -MrlJ't 

0<s<T 



= E [max{Mt, Mt,T + z}\J^t]-E [Mr \ Tt] 

= Mt + E [iMt,T + z- Mt)+ 1 ^t] - E [Mt \ ^t] 

POD 

= E [{MT-t + z- a)+] - / Fr^tix) dx. 

J a 

where a = Mt — Xf. The martingale representation follows from the Clark- 
Ocone formula of Theorem 15.71 □ 

This result extends the martingale representation of the running maximum 
of Brownian motion. 
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Appendix A. Proof of Lemma [5751 

We have that 

sup\\D^^^Fk\\L2{[o,T]xn) < 

k>l 

and 

SUp||L>(2)Ffc||i2([o,T]xRxQ) < OO. 
k>l 

Since L^{[0,T] x Q) and L'^{[0,T] x R x n) are reflexive Hilbert spaces, 
there exist a subsequence (%)j>i, an element a in L^([0,T] x 17) and an 
element (3 in L'^{[0, TjxRxfl) such that D^^^F^. converges to a in the weak 
topology of L^([0,T] x Q) and D^'^'^F^. converges to j3 in the weak topology 
of L2([0, T] x M X S7). Consequently, for any /i G ^^([0, T]), g G L2([0, T] x M) 
and / G L^(S(j^ j,^)(^)), we have that 

{D^'^F,^,h® J^,^_,^){f)) ^^^^^^^^^^^-^ {a,^ 

and 

when j goes to infinity. 

Let F = E[F] + Zn=i ^(n,...,.„)(/^'^'-""^) and F^^ = E[F,J + 

i„) <^(ji,...,i„)(/fc*^' '*"■') be the chaos representations of F and 
Fj;.. By definition, we have that 

(24) dI'^ Fk^ = fl]\t) 

CO n 
n=l (ii,. ..,«„) k=l 

By the linearity of the iterated integrals, the convergence of Ff^.. toward 
F implies that 



n=l (ii,...,i„) 



E E p''"-""' - ft 

n=l 



L2(n) 

2 



in) _ f{il,---,iu) 
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goes to when k tends to infinity. Consequently, it implies that each 

converges to '•••'*") when j goes to infinity. So, using Proposition 14.11 and 
the expression of the derivative in Equation ([241) > g^t that 

n+1 „T 

fc=l-^0 

n+1 .-p 

fc=l-^'^ 

and, as j goes to infinity, this quantity tends to 

n+1 „y 

fc=i-^° 

This holds for any multi-index and functions h and /. Conse- 

quently, 

a(t) = /«(t) 

oo n 
n=l k=l 

and F belongs to D^^^ with D^^^F = a by the unicity of the weak limit. 
Moreover, for any weakly convergent subsequence the limit must be equal 
to D^^^F and this implies the weak convergence of the whole sequence. The 
same argument works to prove that F belongs to D^^^ and that {D'^'^'> Fk)k>i 
converges weakly to D^'^^'F in L'^{X X X P). 

Appendix B. Interpretation oe the directional derivatives 
We consider the product probability space 

which is the product of the canonical space of the Brownian motion W and 
the canonical space of the pure-jump Levy process 

Nt= I I zN{ds,dz) 
Jo Jr 

previously defined in Equation (fTO]l : see Sole et al. [25] for more details on 
this last canonical space. Since L'^{Qw x^n) is isometric to L'^{Qw] L'^{^n)) 
and to L'^{^n'i L'^i^w)) as Hilbert spaces, we will use the theory of the 
Brownian Malliavin derivative and the Poisson random measure Malliavin 
derivative for Hilbert-valued random variables (see [H] and [20]) ■ This is 
possible because both operators are closable. 
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The Brownian MaUiavin derivative for Hilbert-valued random variables 
will be denoted by and the Poisson random measure Malliavin derivative 
for Hilbert-valued random variables by . If we define W = (VFt)tg[o,T] on 
flw X r^jv by 

and N = {Nt)te[o,T] by 

then the process Xt = + aWt + Nt has the same distribution as our 
initial Levy process Xt = fit + aWt + Nt. For notational simplicity, in what 
follows we will write Wt{uj) and Nt{u>') instead of Wt{uj,Lj') and Nt{uJ.,oj') 
respectively. 

We will proceed by induction. If F = f{t) W{dt), then clearly 

Dl'^F = DfF = f{t) and D^F = = 0, 

while if G = g{t, z) N{dt, dz), then 

= D^G = and dJJg = Df^^G = g{t, z). 

Thus, for a fixed n > 1, we assume that -D^^^ and coincide for any 
random variable with chaos expansion of order n. First, let F be of the form 

^= •••®/n®/n+l) = / g{s)fn+\{s)W{ds), 

Jo 

where 

(25) g{s) = J(i,,...,i„) (/i • • • ® /n Is„{s)) • 

To ease the notation, J(ji,...,i„)(/i .../„) will mean - <^ fn)- 

Using the commutativity relationship of Remark 15.31 and the hypothesis of 
induction, we have that 

Z)f V = fn+i{t)git) + U+iis)Dl'^gis) W{ds) 
= fn+i{t)g{t) + r Ui{s)DYg{s) W{ds), 



which is exactly Df^ F, by the classical commutativity relationship of Equa- 
tion 

Secondly, now let F be of the form 

= «^(il,...,j„,2)(/l ® • • • ® /n <X) /n+l) = / / gis-)fn+lis,z)N{ds,dz). 

Jo Jr 

We will use of the integration by parts formula for semimartingales, that is 

[y(i),y(2)]^ = y«y/2)_ f\ii) ^yi2) _ f Y^^J dvP 

Jo Jo 
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if Y^^^ and F^^^ are semimartingales; see Protter [22] for details. If 



(1) 



("2) t ~ 

and Y^ = fn+iis , z) N {ds , dz) , we get that 



F = giT) I I fn+i{s,z)N{ds,dz) 
'0 JR 

I I fn+i{s,z)N{ds,dz)dg{t) 
lo Jo Jm. 



g[-), 1^ I fn+iis,z)N{ds,dz) 



T 



We now consider the two cases where in = 1 and in = 2 separately. We 
have that 



if in = 1; 



1 /o /r Hs-)fn{s, z) N{ds, dz) \iin = 2 



where h{s) = J{i^,...^i^) (/i <8) • • • ® fn-i Is„_i(s)) • If = 1, then 



F = g{T) / / fn+i{t,z)N{dt,dz) 



fn+i{s,y) N{ds,dy) 



h{t)fn{t)W{dt). 



If in = 2, then 



F = g{T) / / /„+i(t,z)7V(dt,dz) 



fn+i{s,y) N{ds,dy) 



T 



h{t-)fn{t,z)N{dt,dz) 



h{t-)fnit, Z)fn+lit, Z) N{dt, dz). 



Note that the last term is an iterated integral of order n (with respect to 
N{dt,dz) for the outermost integral, not N{dt,dz)) since h is an iterated 
integral of order n — 1. So, by the hypothesis of induction, D^^'^ and 
agree for this functional. This is also true for g{T). 
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Consequently, we repeat the previous steps backward with D^^\ If i„ = 1, 
then 



D^F = {dY9{T)) 







h{t)fn{t) 
rT 



D?g{T) 



h{t)Ut) 









fl 

Jo Jm 


D« (g{T) r [ /, 
V Jo JR 




fl 


JO . 


Jo Jm 


[ fn 

Jr 




U 



fn+i{s,y)N{ds,dy) 
fn+i{s,y) N{ds,dy) 
fn+i{r,y) N{dr,dy) 
fn+i{s,y) N{ds,dy) 
fn+i{s,y) N{ds,dy) 
fn+i{r,y) N{dr,dy) 



Df'h{s))fn{s)W{ds) 



Di'h{s)] fn{s)W{ds) 



fn+i{r, y)N{dr,dy) 



h{s)fn{s)W{ds) 



fn+i{r, y)N{dr,dy) 



dg{s] 



and if in = 2, then the same steps are vahd since and D^^^ coincide on 
the extra term. 

The equivalence between D^^^ and follows from the following fact: 
for a fixed n > 1 and a fixed multi-index the linear subspace of 

L^(S(j^^ generated by functions of the form 



(26) 



/i 



fn 



is dense. Indeed, for / G L^(S(jj^ there exists a sequence {fn)n>i, 

whose elements are finite sums of functions as in Equation (f26l) . that con- 
verges to /. We know that D^^^ and are equal for each /„. Since 
and are continuous (see Lemma [531) . they also coincide for /. 

We can apply the same machinery to show that and D^"^^ are the 
same. 
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